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Abstract: We present a novel method for checking the Hurwitz stability of a polytope of 
matrices. First we prove that the polytope matrix is stable if and only if two homogenous 
^ ; polynomials are positive on a simplex, then through a newly proposed method, i.e., the 
weighted difference substitution method, the latter can be checked in finite steps. Exam- 
ples show the efficiency of our method. 
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O ■ 1 Introduction 
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Given a linear time-invariant system x(t) = Ax(t), it is asymptotically stable if the 
system matrix A is Hurwitz stable, i.e., all eigenvalues of A have negative real parts. 
Sometimes, we need to consider a type of system uncertainty in which case the family of 
^ . system matrices forms a polytope, i.e., A is varying in 

m 



A = < <lkA k : ^ Qk = 1, 9fc > for all k 



k=l k=l 



where A\, . . . , A m e ~R nxn are constant matrices. We say that A is robustly Hurwitz stable 
if each matrix in A is Hurwitz stable. The stability of a polytope of matrices cannot be 
derived from the stability of all its edges [1] , which is the case of the stability of a polytope 
of polynomials [2]. In fact, [3] proved that for a polytope of n x n matrices, the stability of 
all In — 4 dimensional faces can guarantee the stability of the polytope, and the number 
2n — 4 is minimal. But checking the stability of 2n — 4 dimensional faces of a polytope is 
also a difficult task. For a matrix polytope with normal vertex matrices, [4] proved that 
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the stability of vertex matrices is necessary and sufficient for the stability of the whole 
polytope. Meanwhile, serval sufficient criterions [5-11] are provided to check stability of 
matrix polytopes. 

Based on the newly proposed results on checking positivity of forms (i.e., homogenous 
polynomials), we present a method for checking the stability of a polytope of matrices in 
this paper, this method is complete, moreover, it only has a power exponential complexity. 

2 Notations 

• JR.: the field of real numbers. 

• Z: the set of all integers. 

• N: the set of all nonnegative integers. 

• M mxn : the space ofmxn real matrices. 

• I n : the identity matrix of order n. 

• det A: the determinant of a square matrix A. 

• Af. the Hurwitz matrix of the polynomial f(x) = a n x n + a n _ix n_1 . . . + ao, it is an 
n x n matrix defined as 
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The successive principal minors of Af are denoted by k = 1, 2, . . . , n. 

• {k\k 2 . . . km): a permutation of {1, 2, . . . ,m}, which changes % to ki, i — 1, . . . , m. 

• m : the set of all ml permutations of {1, 2, . . . ,m}. 

• deg(/): the degree of a polynomial. 

• A T : the transpose of a matrix or vector A. 

• S m : the m — 1 dimensional simplex in W 71 , i.e., 

m 

S m = {Oi, ■ ■ • , x m ) : Xj = 1, Xj > 0, % = 1, . . . , m}. 

i=l 

• [x\: the largest integer not exceeding the number x. 
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3 Main Results 



Suppose 

m 

A = 22 QkA k 
fc=i 

is a matrix in A, denote its characteristic polynomial by 

f A (s) = det(s/ n - A) = s n + a n ^s n ~ l + . . . + a x s + a . (2) 

It is well known that, for all 1 < % < n, (— l)*a n _j equals the sum of all principle minors 
of order i of A, hence form of degree % on qi, q 2 , ■ ■ ■ , q m - Denote by 6,j the (i, j)th 

entry of the Hurwitz matrix of then bij = 0, or bij is a form of degree 2j — i on 

91,92, • • • ,9m- Since 

A fc = ^ ±b lh b 2j2 ■ ■ ■ b kjk , 

and 

bijihh ■ ■ ■ b k j k = 0, 



or 



deg(b lh b 2j2 ■ ■ ■ b kjk ) = ^(2^ - i) = 



k 

i 



k{k + l) 



2 

i=l i=l 

we can see that is a form of degree k{k + l)/2 on q±, q 2 , . . . , g m . 

Theorem 1. There exists a Hurwitz stable matrix in the matrix polytope A, then A is 
Hurwitz stable if and only if 

A n _!(gi, . . . , q m ) > and a (q 1 , q m ) > 0, (q u . . . , q m ) G 5 m . (3) 

Proof. The necessity is directly from Routh-Hurwitz criterion. If A is not Hurwitz stable, 
then by continuity there must exist a matrix A in A which has eigenvalues lying on the 
imaginary axis. Suppose eigenvalues of A are s\, . . . , s n . If some Sj equals zero, then 

a = (-l) n sis 2 • • • s n = 0, 

which contradicts the hypothesis ao > 0. If some Sj and Sj are a pair of conjugate eigen- 
values of A on the imaginary axis, then from Orlando's formula [12], 

A n „ 1 =(-i)^ n ( 8i + 8j )=o, 

l<i<j<n 

which contradicts the hypothesis A n _i > 0. □ 
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Remark: [5] also proved that a polytope matrix A is robustly stable if and only if an 
associated form 7 det A is positive on a simplex, where A is the Kronecker sum of the 
matrix A in A with itself, and 7 is the sign of det A. A is an n 2 x n 2 matrix, hence the 
degree of 7 det A is much larger than A„_i and a in (3). 

A newly proposed method, i.e. the difference substitution method [13,14], can be used 
to check positivity of forms efficiently [15] proved that a form is positive if and only if we 
can get forms with positive coefficients after finite steps of varied forms of substitutions, 
i.e., weighted difference substitutions(WDS). [16] further gave a bound for the number of 
steps required, and pointed out that the WDS method is complete in checking positivity or 
nonnegativity of integral forms. We will introduce this method more detailedly in Section 
4. Based on this method, we have 

Theorem 2. Suppose entries of Ak,k = l,...,m in (1) are all rational, the magnitudes 
of coefficients of A n _i and ao in (3) are bounded by M , then the Hurwitz stability of A 
can be checked by an algorithm with complexity 

O (m m+l n 2m2 {n 2m In M + n 2(m+1) In m + 2(m 2 + mn 2m ) In n) V 



4 Positivity of forms on simplices 



In this section, we will introduce the WDS method for checking positivity of forms [16], 
and analyze the complexity of checking Hurwitz stability of matrix polytopes through this 
method. 

Suppose 9 = (kik 2 . . . k m ) G 6 m , let Pq = (pij) m xm be the permutation matrix corre- 
sponding 9, that is 

1, j = K 

0, j^h ■ 



Given T m G 



where 



/I \ 
i 



T 



m ' 



m 



(4) 



let 



V ... i / 



Ae — PeT„ 



and call it the WDS matrix determined by the permutation 9. The variable substi- 
tution x = Agy corresponding 9 is called a WDS, where x = (xi, £2, • • • , x m ) T , y = 

(yi,V2, ■ ■ -,y m ) T - 

Let /(x) G M[xi, X2, • • • , x m ] be a form, we call 



WDS(J) = |J {/(A,x)} 



(5) 
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the WDS set of /, 

WDSW(/)= |J U {HAe k ---A ei x)} (6) 



t ee m 0iee„ 



the fcth WDS set of / for positive integer k, and set WDS (0) (/) = {/}. 

Let a = a 2 , . . . , a m ) £ N m , \a\ = a± + a 2 + • • • + ct m . For a form of degree d 



f(x!,X 2 , ...,X m )= ^2 



Cry -■- 'y» 



|a|=d 

if all coefficients c Q are nonzero, we say / has complete monomials. 

It is obvious that if there exists fcsN, such that forms in WDS^(/) all have complete 
monomials, and their coefficients are all positive, then / is positive on S m . In fact, the 
reverse is also true, and for integral forms, the upper bound for k can also be estimated. 

Theorem 3 ( [16]). Suppose f £ Z[a?i, x 2 , • • • , x m ] is a form of degree d, and the magnitudes 
of its coefficients are bounded by M, then f is positive on S m , if and only if there exists 
k < C p (M,m,d), such that each form in WDS^ k \f) has complete monomials, and its 
coefficients are all positive, where 



C p (M,m,d) 



In ( 2 dm M dm+l m dm+1+d d^ m+l ^ d+mdm {d + i)( m - 1 )( m + 2 ) 



mm — ln(m — 1) 



(7) 



Remark: The C p (M,m,d) in (7) provides a theoretical upper bound of the number of 
steps of substitutions required to check positivity of an integral form. In practice, numbers 
of steps used are generally much smaller than this bound [15]. 

proof of Theorem 2. A form f(qi,---,q m ) of degree d has at most 

d + m ~ 1 ) <( d + i r -i 
m — 1 / 

monomials thus the number of arithmetic operations of computing WDS(/) is bounded by 

m\(d + l) m -\d + l)^™- 1 ) < m m (d + l) m \ (8) 

Moreover 

C P {M, m,d) = (m(d m In M + d m+1 In m + (m 2 + md m ) In d)) , 

and deg(A n „x) = n{n — l)/2,deg(a ) = n, therefore the complexity of our method for 
checking the robust Hurwitz stability of a polytope of n x n matrices is 

O (m m+1 n 2m2 (n 2m In M + n 2(m+1) In m + 2(m 2 + mn 2m ) In n) V 

□ 
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5 Examples 



First we will illustrate our method through an example from [1]. Suppose A is a 
polytope of following matrices 



A, 



( -1 


-1 

/ -1 



/ -1 



1 







1 \ 

-1 
0.1 







o \ 
1 

0.1 

-1 \ 

-1 

0.1, 



let A = q±Ai + q 2 A 2 + (I3A3, and denote by the characteristic polynomial of A. The 

second successive principle minor of the Hurwitz matrix of is 



63 
25 



99 



243 



144 



qi 



25 50 25 



153 
~25 



qiq2q3 



144 2 243 2 63 , 99 , 171 

q\ q2 + — ^2^3 + rr q2 



q2 2 q3 



25 50 25 25 50 

which is obviously positive on S3. The constant term of Ja{s) is 



qz 



9 o 9 q 23 

«o=T7:gi +tt:<?2 -—q^qs 

10 10 5 

- 3/10 g 2 g 3 2 - 3/10 g 3 V + 



13 2 7 2 



7 

To 



91 ?2 



13 

To 



q2 qz 



19 

To 



(9) 



q3 



do is not positive on S3 since the following form (with a difference of a positive constant 
factor) belongs to WDS^(ao) and its coefficients are all negative: 



3888 qx 2 q 3 - 3888 q 1 2 q 2 - 1568 q 2 q 3 2 



- 6516 qiq 2 q 3 - 1296 gi 3 - 891 q 2 3 

- 2828 q 3 2 qi - 3483 q x q 2 2 - 2223 q 2 2 q 3 - 236 g 3 3 , 

therefore the polytope A is not Hurwitz stable. 

Furthermore, we have checked 900 polytopes of matrices for n = 2, 3, 4 and m — 2, 3, 4, 
i.e., 100 polytopes for each pair (n, m). The vertexes of these polytopes are generated 
following a similar method as was described in [11]: their entries are real numbers with 4 
significant numbers and uniformly distributed in the interval [—1, 1], moreover the maximal 
real parts of their eigenvalues equal —0.0001 (if not so, a shift should be performed). Table 
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Table 1: Time used to check robust Hurwitz stability of 100 polytopes for each pair (n,m) 



n 


m 


number of stable / unstable polytopes 


total time (in seconds) 




2 


67 / 33 


0.125 


2 


3 


29 / 71 


1.578 




4 


11 / 89 


2537.360 




2 


58 / 42 


0.123 


3 


3 


28 / 72 


2.009 




4 


9/91 


665.129 




2 


54 / 46 


0.090 


4 


3 


21 / 79 


2.608 




4 


4/96 


1894.602 



1 shows the time used to check the stability of these polytopes on a computer equipped 
with Intel Core 2 Duo E4500 CPU at 2.2 GHz and 4.0 GB of RAM memory, our program 
have been implemented in the computer algebra system Maple. 

Remark: Generally, the time gets longer with the increase of n or m, but there are 
some extreme examples that very long time may be spent even for small n and m. For 
example, in our experiment, it takes 1129.656 and 1172.234 seconds respectively to check 
the Hurwitz stability of two matrix polytopes generated for (n,m) = (2,4), this makes the 
time corresponding to (n,m) = (2,4) much longer than that corresponding to (n,m) = 
(3,4) or (n,m) = (4,4) in Table 1. 
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